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Abstract: The colliding wind binary (CWB) systems 77 Carinae and WR140 provide unique laboratories for 

4^ X-ray astrophysics. Their wind-wind collisions produce hard X-rays that have been monitored extensively 

^ by several X-ray telescopes, including RXTE. To interpret these RXTE X-ray light curves, we model the 

O wind-wind collision using 3D smoothed particle hydrodynamics (SPH) simulations. Adiabatic simulations that 

^ account for the absorption of X-rays from an assumed point source at the apex of the wind-collision shock 

, cone by the distorted winds can closely match the observed 2-lOkeV RXTE light curves of both 77 Car and 
WR140. This point-source model can also explain the early recovery of 77 Car's X-ray light curve from the 

^ 2009.0 minimum by a factor of 2-4 reduction in the mass loss rate of 77 Car. Our more recent models relax the 

Cn| point-source approximation and account for the spatially extended emission along the wind-wind interaction 

^ shock front. For WR140, the computed X-ray light curve again matches the RXTE observations quite well. 

^ But for 77 Car, a hot, post-periastron bubble leads to an emission level that does not match the extended X-ray 

^ minimum observed by RXTE. Initial results from incorporating radiative cooling and radiatively-driven wind 

^ acceleration via a new anti-gravity approach into the SPH code are also discussed. 

1 Point-Source Emission Model 

77 Carinae and WR140 are long period, highly eccentric colliding wind binaries (CWBs) that provide 
unique laboratories for X-ray astrophysics. Their wind- wind collisions produce hard X-rays that are 
strongly phase dependent; the wide range in stellar separation from the high eccentricity affects the 
X-ray emission, and the severe distortion in the wind geometry from significantly high orbital speeds 
around periastron affects the X-ray absoprtion. Our initial attempts to model the 2-10 keV RXTE 
light curves of both r] Car and WR140 have applied a simple model of point-source emission plus 
line-of-sight wind absorption to 3D, adiabatic, smoothed particle hydrodynamics (SPH) simulations 
of the binary wind- wind interaction (see Okazaki et al. 2008 for details). To match the recent shorter 
minimum of r] Car (Corcoran et al. 2010), the primary mass loss rate is reduced by a factor of 2.5 at 
phase 2.2. Many of the light curve's features are reproduced, including the shorter recent minimum. 
The WR140 light curve matches remarkably well (see Fig. [T] for the light curves and Table [T] for the 
parameters of the SPH simulations). 
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Figure 1: Best-fit RXTE light curves using the point-source emission model for ij Car (left) and 
WR140 (right). The model light curves are normalized to match the average maximum just before the 
first periastron passage (phase ^0.98) and the following apastron (phase 1.5). The free parameters of 
the model- opacity argument of periastron cj, and inclination i- are shown for the relevant models. 



Table 1: Parameters used in the adiabatic SPH simulations of rj Car and WR140. 



Parameter 


rj Car 


WR140 


References 


Primary Mass (M©) 


90« 


50" 


a: Hillieretal. 2001 


Secondary Mass (M©) 


30^ 


196 


b: Marchenko et al. 2003 


Primary Radius (R©) 


90t 


12-^ 


c: Vemer et al. 2005 


Secondary Radius (R0) 


30t 




d: Williams et al. 1990 


Primary Mass Loss Rate (M© yr~^) 


2.5x10-^-'' 


1.2x10-69 


e: Pollock et al. 2005 


Secondary Mass Loss Rate (M© yr~^) 


LOxlO-'^-^ 


3.8 xlO-'^^? 


f: Pittard and Corcoran 2002 


Primary Wind Velocity (km s~^) 


500« 


32009 


g: Zhekov and Skinner 2000 (Model A) 


Secondary Wind Velocity (km s~^) 


3000-'' 


28609 


h: Corcoran 2005 


Orbital Period (d) 


2024^^ 


2899'' 


i: Corcoran et al. 2001 


Semi-major Axis (au) 


15.4* 


16.3" 


j: Dougherty et al. 2005 


Orbital Eccentricity 


0.9* 


0.88" 




Inclination (deg) 


?? 


122-?' 


t : assumed M(M0)/R(R0) =1 


Argument of Periastron (deg) 


?? 


46.7" 





2 Extended Emission Model 

Our more recent efforts to model the RXTE light curves of r] Car and WR140 relax the point-source 
approximation (e.g. Parkin et al. 2009). The emission now originates from spatially extended post- 
shock gas in the wind- wind collision region according to p^K{E,T), where p is the density and 
A(£', T) is the emissivity as a function of energy E and temperature T obtained from the MEKAL 
code (Mewe et al. 1995). The extended wind absorption is now also energy dependent. We then use 
the SPH visualization program SPLASH (Price 2007) to calculate the ray-tracing through the system, 
which generates images in various X-ray bands that combine to make a 2-10 keV X-ray light curve. 
Once again, the WR140 light curve matches well (assuming the opacity is 10 x the opacity of an O 
star wind at solar abundances, an assumption that will be relaxed in future work). The same is not 
true for rj Car, however, where a hot, post-periastron bubble blown into the slow, dense primary wind 
by the much faster companion wind prevents the reproduction of the extended X-ray minimum (see 

Fig.©- 
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apastron and the average maximum just 
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phase 



n Car's 3 keV flux at phase 0.025 
(arbitrary flux units) 

Stars located near 
point of maximum 
X-ray flux, 
primary on right, 
secondary on left 

^hot, post-periastron 
bubble 



Figure 2: WR140 RXTE light curve using the extended emission model (left) and the post-periastron 
bubble of t\ Car (right). 



Radiative cooling, via the Exact Integration Scheme (Townsend 2009), and the radiative driving of 
the stellar winds, via an anti-gravity approximation (an assumed outward acceleration that competes 
with gravity instead of the full Castor- Abbott-Klein calculation), have been implemented. The accel- 
eration of the secondary wind towards the primary star drastically decreases around periastron due 
to radiative inhibition, where the radiation from one star provides enough force on it's companion's 
wind to inhibit it's acceleration (Stevens & Pollock 1994). This leads to the secondary's wind collid- 
ing with the primary's wind at a much lower speed directly between the stars, so the high temperature 
shock cone of t\ Car collapses, as was proposed by Parkin et al. (2009). However, there is still a high 
temperature shock on the backside of the secondary, the side opposite the primary star. The high 
orbital speed at periastron essentially causes the secondary to become buried in the primary's wind, 
so the secondary wind on the backside will collide with primary wind material. Since the acceleration 
of the backside's secondary wind is not affected by radiative inhibition, this portion of the wind still 
collides at a high speed (although the collision is co-moving instead of head-on) to create a high tem- 
perature, post-periastron bubble. The hard X-rays produced in this post-periastron bubble prevent the 
reproduction of the minimum of the RXTE light curve with the extended emission model. Outside 
the minimum, the model matches fairly well (see Fig. [3]). 

3 Conclusions & Future Work 

The results of using 3D SPH simulations to model the wind- wind interaction of the CWBs 77 Carinae 
and WR140 have been presented, as well as our attempts to model these systems' RXTE light curves. 
Although the point-source emission model produces remarkably good results for both systems when 
applied to the adiabatic SPH simulations, this model is clearly an approximation that needs to be 
relaxed. Upon doing this by implementing an extended emission model, the resultant X-ray light 
curve for WR140 matches pretty well, while t] Car's light curve does not match well during the X-ray 
minimum. A hot, post-periastron bubble forms in Car from the fast companion wind blowing into 
the slow primary wind just after periastron, and this large emitting region creates a significant digress 
from the point- source result. On the other hand, the similar wind speeds in WR140 prevent a post- 
periastron bubble from forming, so the emission is approximately point-like throughout the orbit, and 
hence in good agreement with the point-source emission model results. 

To improve the match of the extended emission model during r] Car's RXTE X-ray minimum, 
radiative cooling and radiative driving of the winds via an anti-gravity approach were added to the 
SPH code (previous SPH simulations were adiabatic and had the winds leave their star at terminal 




Figure 3: Density and temperature snapshots of an SPH simulation of r] Car that includes radiative 
forces and radiative cooling (left) and the resulting RXTE model light curve using the extended emis- 
sion model (right). In the SPH snapshots, r] Car is the white dot, while the companion is shown by 
the black dot in the left density panel and then orbits counterclockwise. The snapshots show that hot 
material exists between the stars before periastron (left panels), vanishes at periastron (center panels) 
from radiative forces and radiative cooling, and begins to reappear after periastron (right panels). The 
right panels also shows the hot, post-periastron bubble that forms below both stars, which prevents 
the matching of the RXTE light curve during the X-ray minimum. 

speed). While this leads to the collapse of the high temperature shock between the stars at periastron 
due to radiative inhibition, it does not affect the hot gas in the remainder of the post-periastron bubble. 
Therefore, there is still excess X-ray emission at periastron compared to the RXTE observations. 

At the time of this conference, the results from the SPH simulations with radiative cooling and 
radiative driving are relatively new, so our immediate future work will be to further explore and under- 
stand these results. Of particular interest is the post-periastron bubble, which appears to be preventing 
the matching of the RXTE light curve during the X-ray minimum. Furthering the intrigue about this 
bubble are the results of Parkin et al. (these proceedings) where a 3D grid-based hydrodynamics code 
was used to model 77 Car, and although the density and temperature structures compare very well from 
the grid-based code to the SPH code, and they too have excess emission during the X-ray minimum, 
the authors do not claim to see a post-periastron bubble. Other future work will include modeling 
WR140 with the radiative cooling and radiative driving improvements to the SPH code, modeling 
other CWBs such as HD5980, and improving the resolution of the extended emission model to begin 
modeling the X-ray spectra of CWBs. 
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